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(g) Asservlssenr)ent de phase, en bande de base. 

(57) Le proc^d6 consiste d extraire (5) la phase de 
chacun des signaux complexes en supprimant 
las sauts de phase inconnus, d rechercher (6) 
dans la suite des phases obtenues une premiere 
plage de phase de longueur d^termin^e L pour 
laquelle la variation de phase est la plus lin^aire 
possible, k d6terminer (7) une deuxi6me plage 
de phases de longueur L2 sup^rieure ^ L conres- 
pondant d une approximation lin^aire plus pr^ 
cise. ci effectuer (8) it^rativement des 
regressions Iin6aires d I'aide tfune fenetre glis- 
sante de longueur 6gale ci la deuxieme plage de 
phase L2 pour supprimer les sauts de phase 
r^siduels pour effectuer (9) une regression qua- 
dratique sur les phases restantes et pour calcu- 
ler une estimation globate des phases de 
reference des signaux de cheque bloc, et k 
effectuer (10) une comparaison de I'enreur qua- 
dratique moyenne entre les phases de reference 
estimSes et les phases brutes obtenues par 
regression lineaire. 

Application : Transmission de donnSes par 
satellite defilant. 
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La pr6sente Invention concerne un proc6d6 et un dispositif de d^niodulation coh^rente par blocs de si- 
gnaux complexes niodul^s en phase. 

Elle s'appllque aux transmissions de donntes num^riques par vole radlo^lectrique avecfalble rapport si- 
gnal ^ bruit et dScaiage en frequence variable dans le temps et notamment aux liaisons radlodlectriques entre 
stations au sol via un satellite d^filant. 

Dans les syst^mes de demodulation coh^rente de modulations ^ 2 ou 4 6tats de phase, il est connu, pour 
am^liorer les performances de transmission de messages de taille minimale en prteence d'un fort niveau de 
bruit et d'un fort d^calage en frequence non constant (rampe Doppler), d'utitiser un codage de type convolutif 
consistant ^ 6mettre plus de symboles que ce quil est strictement necessaire. Les codes utilises le plus sou- 
vent sont des codes de redondance R = 2 ce qui conduit k un d^bit utile 6gal k la moiti6 du d^bit transmis, un 
bit d'information se traduisant par 2 bits transmis. 

En throne, le rapport signal/bruit ^nergetique par bit utile ndcessaire pour obtenir un taux d'erreur de 1 0-^ 
est de I'ordre de 6 dB avec une marge de realisation de 2 dB. Dans ces conditions, le rapport signal/bruit en 
entr6e d'une demodulation de reception est de 6 dB - 10 logioR soit 3 dB, ce qui est tr^s faible. Cela impose, 
dans les systemes dassiques de faire pr^ceder tout message d'une sequence de synchronisation initiale suf- 
f Isamment longue pour que, m§me dans des conditions de propagation tres defavorables, la demodulation 
se fasse dans de bonnes conditions. 

Parmi les types de demodulation connus deux apparaissent utilisables, la demodulation dite "au fil de 
Teau'* et la demodulation par bloc. 

Dans la demodulation "au fil de Teau" des echantillons du signal sont fournis en continu ^ Tentree du sys- 
teme de demodulation. 

Le systeme de demodulation cc ,:porte classiquement, d'une part, une boucle d verrouillage de phase, 
destinee k estimer en continu la phase qu'aurait le signal en I'absence de modulation, pour donner une phase 
de reference, et d'autre part, un demodulateur, fournissant en sortie soit directement les valeurs de Tlnforma- 
tion decodee soit le signal complexe redresse. Dans ce systeme, Torgane le plus Important est la boucle d 
verrouillage de phase, dont les caracteristiques sont le resultat d'un compromis entre la vitesse d'estlmatlon, 
dont depend le temps d'accrochage, et la "gigue" qui conditionne le taux d'erreur residue!, ou dans le meilleur 
des cas, la degradation des performances par rapport au cas ou la phase du signal serait parfaitement connue. 

Dans la demodulation par blocs, des echantillons complexes du signal sont d'abord stockes dans une me- 
molre puis, quand lis sont en nombre suffisant, ceux-ci sont traites en bloc en une ou plusieurs passes. La 
construction des demodulateurs par blocs est recente car elle a ete rendue possible par les progrds techno- 
logiques en matiere de memoires et de microprocesseurs de traitement de signal, seuls capables d'effectuer 
une demodulation complete sur les echantillons stockes en memoire. Dans les systemes S demodulation par 
blocs la demodulation a lieu en deux operations successives, une premiere operation qui consiste e estimer 
la phase de reference du signal par mise en oeuvre d'algorithmes divers travaillant sur la totallte des echan- 
tillons stockes en memoire et une deuxieme operation consistant ^ corriger la phase des echantillons et decider 
la suite du traitement La suite des traitements consiste par exemple ^ effectuer un decodage au maximum 
de vraisemblance, le plus souvent en faisant usage de I'algorithme de Viterbi, lorsque le signal comporte un 
code convolutif. La meme structure de materiel est utillsee pour effectuer d la fols la demodulation et le de- 
codage. 

Une description d'un tel systeme, non compris le decodage correcteur, et la comparaison de ses perfor- 
mances avec celles d'un systeme plus classique sont decrits dans I'artlcle intitule "Block Demodulation for 
Short Radio Packets'* de J. Namiki, Electronics and Communications in Japan, J67-B, N'^l (Janv. 1 g84) pp.54- 



Comme indique dans I'article "Timing Extraction and Carrier Estimation for PSK Signal Block Demodula- 
tion System" de Y. Higashida, M. Hagiwara, N. Nakagawa, Electronics and Communications in Japan, J71-B, 
N''4 (Avr. 1988) pp.540-546, les systemes de demodulation orecedemment decrits ont d'enormes diff icultes 
pour fonctionner a de tres faible rapport signal/bruit et une evaluation des performances de ces systemes mon- 
tre clairement que le rapport signal/bruit minimal de fo notion n erne nt est plus proche de 10 dB que de 3 dB. 
Qui plus est, les cas ou le decalage Doppler varie au cours du temps ne sont pas traites. 

Le but de Tinvention est de pallier les Inconvenients precites. 

A cet effet. Tinvention a pour objet un precede de demodulation coherente par bloc de suites de signaux 
complexes Sn modules en phase sulvant des etats de phase determines caracterise en ce qu'il consiste a ex- 
traire la phase de chacun des signaux complexes en supprimant les sauts de phase inconnus, e rechercher 
dans la suite des phases obtenues une premiere plage de phase de longueur determlnee L pour laquelle la 
variatton de phase est la plus lineaire possible, e determiner une deuxieme plage de phases de longueur L2 
superieure e L correspondent d une approximation lineaire plus precise, e effectuer iteratlvement des regres- 
sk)ns llneaires ^ I'alde d'une fen§tre gllssante de longueur egale k la deuxieme plage de phase L2 pour sup- 
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primer les sauts de phase residuals pour effectuer une regression quadratique sur les phases restantes et 
pour calculer une estimation globale des phases de r6f6rence des signaux de chaque bloc, et d effectuer une 
comparaison de Terreur quadratique moyenne entre ies phases de reference estim^es et les phases brutes 
obtenues par regression Iin6aire. 
5 LInvention a pour principal avantage qu'elle permet, la recherche de la phase de reference avec une faible 

charge de calcul. Elle permet 6gaiement de rendre minimum la degradation des performances pour les tr6s 
faibles rapports signal/bruit, en penmettant une qualification de la quality du signal regu et d'aff iner la syn- 
chronisation. 

D'autres caract6ristiques et avantages de I'invention apparattront ci-apres d I'aide de la description et des 
10 dessins annexes qui represented : 

- La figure 1 un dispositif de modulation selon rinvention. 

- La figure 2 les etapes principales du procede selon I'inventbn mises sous la forme d'un organigramme. 

- Les figures 3A et SB des diagrammes comparatifs montrant les effets de la correction de phase selon 
rinvention. 

15 - Les figures 4A et 4B des diagrammes de phases conig6es et d*erreur minimum pour un rapport signal/bruit 

de 3 dS. 

- La figure 5 un exemple de droite de regression llneaire traces sur une plage etendue relatlvement d une 
plage initiate. 

- La figure 6 une illustration du processus d'estimation finale des phases de reference par prolongement 
20 d droite et e gauche de la droite de regression. 

- La figure 7 une representation de phases corrigees initlales relatlvement d une regression quadratique 
des phases corrigees pour un rapport S/B = 3dB. 

- La figure 8 une representation de phases corrigees relatlvement e une regression quadratique des pha- 
ses corrigees pour un rapport S/B de 0 dB. 

25 - Les figures 9A, 9B, 9C, 9D des diagrammes de demodulation de signaux complexes e deux etats de 

phase pour des rapports S/B de 3, 0, 6 et 10 dB. 
Le dispositif de demodulation selon invention qui est d6crit ci-apres a une structure demodulateur par 
blocs. II comporte un dispositif 1 de passage du signal S e demoduler en bande de base, un dispositif de me- 
morisation 2, un dispositif d'estimation de phase de reference 3 et un dispositif de demodulation et de deco- 
30 dage 4. Le dispositif 1 extrait de fa^on connue sur une frequence intermediaire Fl les composantes reelles et 
imaginaires du signal complexe S e demoduler. .Ce signal est forme par des suites de N echantillons Sn 
complexe de la forme : 



35 



Sn = a exp j Go + 2 ji Af n T + 2 n dAf 



expj(<Pn) + Zn,n = 1...N 



ou 

40 a = amplitude 

Oo = phase tnitlale 

Af = ecart en frequence initial 

dAf = variation de recart en frequence 

(p„ = phase porteuse de la modulation (M=2 ou 4 possibilites) 
45 T = duree du symbole et intervalle d'echantillonnage 

Zn = bruit gaussien complexe de puissance 2a^ 

a2 

rapport signal/ bruit S/B = — r 

Le dispositif de memorisation 2 est couple en sortie du dispositif 1 et stocke "au vol" chaque bloc de N 

50 echantillons avant tout traitemenL 

Le dispositif 3 estime la phase de reference du signal qui ne depend pas de rinformation transmise. k 
partir des echantillons stockes dans le dispositif 2. Le dispositif 4 est couple d'une part, en sortie du dispositif 
de memorisation 2 et d'autre part, en sortie du dispositif 3 d'estimation de phase pour reconstituer des echan- 
tillons complexes redresses et reconstituer rinformation utile transmise e partir des echantillons complexes 

55 stockes dans la memoire 2 et des phases de reference fournies par le dispositif 3. Le dispositif 3 d'estimation 
de phase opere suivant les etapes 5 d 10 du precede represente e la figure 2. L'etape 5 consiste e calculer la 
phase du signal complexe S applique d I'entree du dispositif 1 et conrige les phases en supprimant des sauts 
de phase inconnus dus e rinformation qui le module. 
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L6tape 6 conslste d rechercher dans la s6rie des phases d6livr6es d I'^tape 5 un plage de phase de lon- 
gueur Lou la variation de phase parait etre la plus lin^aire possible. 

Al'^tape 7 une plage de phases de largeur L2 supSrieure d L est recherch6e correspondant d une approxi- 
mation lin^aire plus precise. 

L6tape 8 effectue it6rativement des regressions lin§aires d I'aide d'une fenetre glissante de longueur L2 
pour supprimer tous les sauts de phase r^siduels. 

L6tape 9 effectue une regression quadratique sur les phases obtenues d retape 8 pour calculer une es- 
timation globale des phases de reference sur chaque bloc de N 6chantillons trait6s. Enfin retape 10 effectue 
une comparaison de f erreur quadratique moyenne entre les phases de reference estim^es d I'etape 9 et les 
phases brutes obtenues ^ retape 8, d une valeur de seuil d6termin6e. 

La correction des phases qui est effectu6e d I'^tape 5 vise d lafois d supprimer influence de la modulation, 
qui se traduit par une incertitude de n sur la phase de r6f6rence du signal et ^ supprimer I'effet des Evaluations 
de phase qui ont toujours lieu d 2it pr6s. Elle s'effectue elle-meme en trois stapes de traitement. La premiere 
etape consiste d corriger les phases des echantillons en les ramenant par example dans I'intervalle 

[It Til 

pour une modulation biphase 0 - n. Dans ce cas, si P„ = arg(Sn) est la phase comprise entre -n et +n du n« 
echantillon de signal, la correction a lieu en transformant Pp pour n = 1 ... d N par les iterations suivantes : 

- faire Qn=Pn 

- si Qn > I alors faire Qn = On - 

- si Qn<-^alors faire Qn = Qn + Jc 

La deuxieme etape consiste k calculer les differences de phase successives et e les ramener elles aussi 

entre -^et + ^. 
2 2 

Ualgorithme correspondant consiste ^ poser 

Dn = Qn - Qn . 1 et pour n = 2 ... N 

e faire D„ = □„ - « si D„ < - | ou ^ faire D„ = D„ + jc si D„ < - 1. 

La troisieme etape vise ^ reconstituer une suite de valeurs de phases On representative (d une incertitude 
de n pres) de la phase de reference du signal. 
L'algorithme correspondant consiste a poser : 

a>i = Qi et ^1 faire On = On - 1 + Qn 

pour n = 2 ... N. 

Des valeurs respectives de Pn et On, pour un bon rapport signal/bruit de 20 dB et pourun mauvais rapport 
signal/bruit de 3 dB sont tracees respectivement sur les graphes des figures 3A-3B ou les phases orlginales 
en radians de N = 100 echantillons sont tracees en traits pointilies et les phases conrigees apparaissent en 
trait plein. 

Comme cela apparart sur la figure 38, les phases apparaissent beaucoup moins bien corrigees lorsque 
le rapport signal/bruit est de 3 dB, car les phases corrigees qui devraient en principe etre quasi continues, pre- 
sentent des sauts de phase residuals de n radians. Cependant. la figure 3B permet d'observer que les phases 
corrigees qui presentent des sauts de phase de n comportent aussi des plages sans sauts. De plus, si le de- 
calage en frequence ne varie pas trop vite, la phase varie (presque) lineairement sur une telle plage. Uint6ret 
de cette constatation a de importance, car elle permet lorsqu'une plage sans sauts est reperee, de corriger 
les sauts de phase en la prolongeant de part et d'autre. Cette particularite est mise d profit par I'invention e 
retape 6 en recherchant dans la suite des phases corrigees, une plage de L valeurs successives approxima- 
tivement reparties sur une droite et pouvant etre approximees par une equation de la forme : 

0^ + j«Ai+ B pouri = 1,..., L 
no etant le numero de rechantillon precedent immediatement le debut du bloc. 

Par regression lineaire le minimum de I'erreur quadratique moyenne E donnee par 
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E(no) = ^i:Ko+i-(Ai + B))' 



1=1 

est calculi. 

Les valeurs A et B correspondant k ce minimum sont donates par les relations : 

= Pi ^0 ~ Pq ^1 
A 

B ~ 

avec 

1=1 ^ i=1 

A = So 82 - 32 

L L 

Po=Z^no+i.Pl=Z'^no+i 
i=1 i=1 

La valeur mininr)ale de Terreur est alors : 



Emin(no) = ^[z<.i-APl-BPo 



II suff It alors de calculer Emin(no) pour no = O...N-L et de rep^rer la valeur de Hq correspondant au minimum. 
La valeur de no d^signe la position du d^but de la plage oCi L ^chantillons successlfs de phase corrigSe sont 
le mieux approch^s par une droite. Pour moins d6penser en puissance de calcul les diff6rentes sommes faisant 
intervenir les phases peuvent aussi etre Svalu6es de proche en proche. La valeur de L doit alors etre soi- 
gneusement choisie d'une part, suff isamment faible pour que, meme avec un faible rapport signal/bruit, il soit 
quasiment certain de trouver une plage de phases corrig^s quasi lineaires, sans saut de phase et d'autre 
part, suff isamment ^lev^e pour que la regression lin^aire donne un r^sultat suffisammentf iable. Pourun rap- 

^ port signal/bruit minimum de 3 dB, une valeur de L de Tordre d'une dizaine d'6chantillons donne des resultats 
acceptables, comme le montrent les graphes de phase corrig§e en radian et d'erreur quadrattque moyenne 
des figures 4A et 4B ou il apparaTt que sur N = 100 ^chantillons, i'erreur est minimale avec un point de depart 
no (offset) de 56 echantillons ^ partir de I'origine. Dans ces conditions c'est-d-partir de cette portion du signal 
qu*est cherchee Testimation de la phase de reference. 

^ La largeur de la plage L est d^finie pour avoir une probability sufftsante de trouver una plage correcte. 

Cependant, comme les valeurs des coefficients de regression lineaire A et B correspondant au minimum 
d'erreur sont entach^s eux aussi d'erreurs, ceux-ci sont estim^s avec plus de precision en cherchant d eiargir 
la plage L sur une longueur L2. La longueur L2 est fonction de la valeur de la rampe Doppler maximum admis- 
sible. Elle est choisie de telle sorte que I'erreur de linearite reste negligeable. inferieure par exemple d la nnoitie 

^ de recart-type de I'erreur due au bruit 

Cependant comme reiargissement de la plage initiate a pour les faibles rapports signal/bruit toutes les 
chances de conduire d prendre en compte des phases corrigees pr^sentant des sauts de phase nmjittples de 
n. il est n^cessaire pour s'accommoder de cette situation d'effectuer les operations suivantes. 

Tout d'abord. la plage en cours d'eiargissement est suppos^e avoir une largeur £ superieure ou egale d L 

^ et inferieure e une largeur L2. Pour reiarglr d'une unite, les operations consistent e calculer une valeur extra- 
poiee par approximation lineaire d partir des valeurs courantes de A et B. e modifier ensuite la phase corri- 
gee On, d prendre en compte en lui ajoutant ou en lui retranchant ntds le nombre de fois necessaire pour que 
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la difference I Op - Oe I soit minimale. La nouvelle valeur de On ainsi obtenue est alors utiliste pour remettre 
^ jour les coefficients de regression lin6aire Aet B. 

Le processus d'6largissement peut commencer par exemple par un prolongement de la plage vers la gau- 
che pour ensuite la prolonger vers la droile. L'algorithme k utiliser peut alors §tre le suivant : 

Prolongement 6 gauche : 

INITIALISATION : 

Remettre les quantlt6s Sq, Si, S2, Pq, Pi. A et B aux valeurs trouv^es ^ la position n©. 
Fairs n = no (indice du premier 6chantillon k considdrer d gauche) 

Faire i = 0 (i vaudra 0, -1 . -2 -(L2-2)/2) 

Fa ire 

ng = MAX(1. M IN(no - (L2 - L) / 2 + 1. N - (L2 - 1))) 

(ce sera le dernier Indice utilise). 
ITERATION : 



Si n < ng alors Fin (prolongation termin6e) 
= = A i + B (extrapolation) 

Ramener le plus pr6s possible de C>e par ajout(s) ou retrait(s) de n 
(mise d jour des coefficients de regression) 

So = So + 1, Si + I. S2 + i^ 
A = S0S2 - S2 

PO = Po + ^n, Pi = Pl = i <»„ 

A = ^1 ~ Pq Si 

A 

Pi Si 



B = ^0 ^2 



e : n = n - 1 
i = i- 1 
Aller en a 

Prolongement e droite 

INITIALISATION : 



40 



45 



50 



55 



Remettre les quantit6s So, Si, S2, P©, Pi. A et B aux valeurs trouvees e la position no- 
Faire n = n© + 1 (indice du premier echantillon a constderer e droite) 
Faire i = L + 1 (i vaudra L + 1, L + 2 L +{L2-L)/2) 

ITERATION : 



Si n > ng + L2 - 1 alors Fin (prolongation terminee) 
<I>e = A i + B (extrapolation) 

Ramener le plus pres possible de Oe par ajout(s) ou retrait(s) de n 
(mise ^ jour des coefficients de regression) 

So = So + 1. Si + i, S2 + |2 

A = So S2 - SJ 

Po = Po + On, Pi = Pi + i On 
Pi So - Pq Si 



A = 



B = 



P0S2 - Pi Si 



n = n + 1 
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i = i+ 1 
Aller en a 

En fin d'ex^cution de ralgorithme. les coefficients de regression lin§alre A et B valables pour une plage 
de largeur L2 contenant la plage initlale de largeur L sont disponibles. 

La figure 5 montre les droites de regression lin^aire obtenues dans le meme referentiel que celul des fi- 
gures 3A et SB correspondant aux deux plages de I'exemple prScMent (la plage initiale L 6tant rep^rte par 
deux traits en tirets paralieies). 

L'estimation finale de la phase de reference a lieu par I'execution des stapes 8 et 9. Comme la longueur 
L2 a §t6 choisie la plus grande possible, le reste des estimations de phase est execute ^ Taide d'une fen§tre 
giissante de largeur L2 pour corriger les sauts de phase. La fenetre trouv6e de largeur L2 trouv6e d retape 7 
Incorpore les 6chantillons de rang ng 6 ng + L2- 1 si L2 est inf6rieurou 6gal au nombre N d*6chantillons. L*6tape 
8 consiste k faire glisser la fenetre k gauche puis k droite. Ce glissement s'ef fectue suivant trois modes op^- 
ratoires consistent, k calculer une valeur nominale extrapol^e <I>o de la phase, k effectuer une correction de 
chaque phase <!>„ et k effectuer une mise k jour de la regression Iln6aire. 

Pour le glissement k gauche ralgorithme k ex^cuter est le suivant : 



INITIAUSATION : 

Remettre quantit6s So, Si. S2, Pq, Pi, Aet B aux valeurs trouvtes k la fin de determination initiale de la 
20 fenetre de largeur L2. 

Leurs valeurs sont : 

= Pi Sq ~ Pq ~ Sj 



25 



3 ^ P0S2 - Pi Si 



A 

"a" 



avec 



30 



_^i_ L2(L2-M) o^ ^j2 >2('-2'^l)(2L2-H) 



i=1 



i=»1 

A = S0S2 - S2 
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ng+L2-1 



Po= Z ^i.Pi= Z ('-ng+l)<D, 



i=n„ 



ig i=ng 

Paire n = ng - 1 (indice du premier echantillon k considerer k gauche) 
ITERATION : 



45 



50 



Si n < 1 alors Fin (plus d'echantillons k gauche) 
Oe = B (extrapolation, pour i = 0) 

Ramener On le plus prks possible de Oe par ajout(s) ou retrait(s) de n 

(mise k jour des coefficients de regression) 

A - Pj Sq - Pq Si 

^ " A 

g ^ P0S2 - PiS, 



avec 



55 



n+L2-1 



n+Lo-1 



Po= S <J^.Pl= Z (i-n + l)<I>i 



i=n 



i=n 
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e : n = n - 1 
Aller en a 

Pour le glissement k droite ralgorithme d ex^cuter est le suivant : 
5 INITIALISATION : 

Remettre quantit^s So, S,, S2, Po. Pi, Aet B aux valeurs trouv6es ^ la fin de determination initiale de la 
fendtre de largeur L2. 

Faire n = ng + L2 (indice du premier tehantillon k consid^rer k droite) 

10 

ITERATION : 

a : Si n > N alors Fin (plus d'6chantillons d droite) 
b : a>e = A (L2 + 1) + B (extrapolation, pour I = L2 + 1) 
15 c : Ramener On le plus pr^s possible de Oe par ajout(s) ou retrait(s) de «rds. 
d : (mise d jour des coefficients de regression) 

A = Pj Sq - Pq Si 
A 

R - Pq S2 - Pj Si 

20 B 

avec 



25 



n+L2-1 n+L2-1 
Po= Z ^i'P^= Z (i-n + l)Oj 
i=n iasn 



30 



e : n = n - 1 
Aller en a 

Cependant pour ne pas dSpenser en puissance de calcul, les calculs precedents peuvent encore etre exe- 
cutes par un ajustement des dif ferentes variables comme def ini ci-apres. 
L'algorithme pour le glissement e gauche est alors le suivant : 



35 



40 



INITIALISATION : 

Remettre quantites Sq, S^, S2, Po, Pi, A et B aux valeurs trouvees S la fin de determination initiale de la 
fenetre de largeur L2. 

Faire n = ng -1 (indice du premier echantillon e considerer ^ gauche) 
Faire i = n - no 

ITERATION 



45 



50 



55 



Si n < 1 alors Fin (plus d'echantillons d gauche) 
a>o = A i + B (extrapolation) 

Ramener le plus pres possible de O^e par ajout(s) ou retratt(s) de nrds 
(mise a jour des coefficients de regression) 

Si = Si + i - (i + L2) = Si - L2 
S2 = S2 + i2 _ (i + L2)2 = S2 - L2 (L2 + 2i) 
A = S2S0 - S2 

Po = Po + On - On ^ L,. Pi = P^ + i On - (1 + L2) 0„ h- ^ 
^ _ Pi Sq - Pq Si 



Q _ Pq S2 - Pi Si 



A 

"a" 



n = n - 1 

i = i- 1 

L'algorithme pour le glissement ^ droite est le suivant : 
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INITIAUSATION : 

Remettre les quantites So, St, S2, Po. Pi. A et B aux valeurs trouv6es k la fin de determination initiale de 
la fen^tre de largeur L2. 

Fa ire n = ng -(> L2 (indice du premier ^hantillon ^ consid^rer k droite) 
Faire i = n - n© 

ITERATION 
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IS 



20 



25 



30 



a : Si n > N alors Fin (plus d'^chantlllons d gauche) 
b : a»e = A i -I- B (extrapolation) 

c : Ramener On plus pr^s possible de Oe par ajout(s) ou retrait(s) de n 
d : (mise d jour des coefficients de regression) 

Si = Si + i - (i + L2) = Si - L2 
S2 = S2 + P - (i + L2)2 = S2 - L2 (L2 + 2i) 
A = S2S0 - S2 

Po = Po + <l>n - ♦ L,. Pi = Pl + i^n " (* + Lj) ♦n * L, 

A - Pi So - Po St 
^ ■ A 

a _ Pq S2 - Pi St 

® A 

e : n = n- 1 
i = i-1 
Aller en a 

Les approximations lin§aires locales de largeur L2 ainsi obtenues permettent d'obtenirun jeu de N phases 
corrigees On ou il n'y a plus de sauts de phase comme cela est represents d la figure 6. 

Le dernier traitementconsiste, compte tenu du fait que le decatage Doppler varie dans te temps, e effectuer 
une regression quadratique sur les phases corrigees, oe qui suppose que le decalage Doppler varie quasi- 
lineairement. 

Ce traitement consiste d approximer les phases par : 

On«An2 + Bn + C. n = 1...N 
et d minimiser Terreur quadratique moyenne suivante : 



N 



E=^Z(*n-(An2.-Br.-hC)) 
^ n=1 



40 



45 



La solution est donnee par : 

P^ Sl - Si S3) + Pi(So S3 - Si S2) P2(S^ - S0S2) 

A 



A = 



_ Po(Si S4 - S2S3) + P^(Sl - S0S4) + P2(SoS3 - St S2) 
B = - 

^ _ PoCSj - S2S4) -I- Pt(St S4 - S2S3) -i- P2(Si - St S3) 



avec 



50 



So = N.Si = X' = 



^. N(N + 1) ^ N(N + 1)(2N+1) 

- -z .S2 = 2-' = H 



i=1 



i=1 



55 



S3 = = 
i=1 



N2(N + lf N(N + 1)(2N + 1)(3N2 + 3N-i) 

■ " -1 — ■S4 = Z' =- — — 



i=1 
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A = S| - 2S1S2S3 + SoS| + S?S4 - S0S2S4 
N N N 

n=1 n=1 n=1 



et la valeur minimale obtenue est : 

r N 



E --I 



Z^n-AP2-BPi-BPo 



Les phases corrig^es initiales, et Tapproximation parabolique finale obtenues parle dernier traitementsont 
representees figure 7. 

Le coefficient B est §gal S rincr^mentation de phase estimte pendant la dur^e d'un symbole, pendant le 
symbole (f ictif) d'indice 0, soit 2n Af T. 

Sa variance est de I'ordre de 192 (a/a)2 /N^ (en radians), a 6tant I'amplitude du signal et a r6cart-type 
du bruit. Le rapport S/B est donnd par : S/B = /{2g^). 

L'^art-type est approximativement ^gal d 



a^f = (1 / 2(2 71 T)) (96 / (n^ S / b)) 



1/2 



A titre d'exemple, pour 1 Kb/s, S/B = 3 dB et N = 2*100 bits transmis (100 bits utiles), a^f estde I'ordre de 0.4 
Hz. 

Cet exennple montre que, lorsqu'aucune erreur n'a commise, la valeur minimale de I'en-eur quadratique 
moyenne estfaible. Pour un rapport S/B minimum de 3 dB, elle est en g6n6ral inf6rieure k 0.2 radian^. 

Par centre, lorsqu'une erreur est commise par une correction abusive de Ttrds, cette valeur minimale est 
proche de ou sup6rieure d la valeur thdorique 7r2/i2=0.82 radian^, qui correspond d une erreur uniformdment 
r^partie entre -7c/2 et +n/2. 

Cette constatatlon penmet de s6lectionner les blocs regus en fonctron de Emin- Si cette quantity est infe- 
rieure d 0.3 (par exemple), le bloc est declare acceptable, sinon it est rejet6. Naturellement ce traitement peut 
tou jours etre recommence avec une valeur de L plus faible (par exemple moitie moindre) af in d'augmenter les 
chances de trouver une plage de longueur L sans transition. La selection penmet aussi d'^liminer les messages 
ayant un trop mauvais rapport signal/bruit, comme le montre la figure 8 oCi S/B vaut 0 dB avec une valeur de 
Emm sup6rieure d 0.3 (le signal d'entr^e §tant le meme que precedemment, mais ou la puissance du bruit a 
ete multlpliee par 2). 

La phase finale de demodulation consiste d soustraire la phase de reference du signal regu et done a cal- 
culer : 



= Sn exp j(-(A + B n + c)). n = 



Une premiere constellation d'^chantillons complexes ainsi obtenues par demodulation sur deux voies en 
quadratures est repr6sent6s k la figure 9A dans le cas d'un rapport signal ^ bruit SNR = 3 dB. 11 n'y a qu'un 
bit errone pour un ^cart moyen de 0,248943.. 

D'autres constellations obtenues avec des rapports S/B de 0 dB (moins bon), puis 6 et 10 dB (done meil- 
^ leurs) sont montr^es respectivement aux figures 9B, 9C et 9D. 

Dans les calculs precedents, les variances des coefficients de regression dependent de la fagon dont la 
phase <I>n est approximee. 

Dans I'approximation quadratique ou est de la forme 0„ ^A^n + Ban + C3 + £„ les variances des coef- 
ficients A3, B3 et C3 sont donnees par les relations : 

V(Ai^18g 
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V(B) 192 

et 

5 c2 N 

<Te repr6sente r6cart type et est voisin de — 1/(2S / B)^'2 

Lorsque Tapproximation de On suit une loi de regression lin^aire de la forme : 

On = Ajn + B2 + en 
10 i^es variances sur les coefficients A2 et B2 sont donn^es par ies relations : 

o2 N3 

et 

IS a2 N 

Dans fapproximation de On par une constante de la forme : 

On = A, + 8b 

on obtient : 

V(A) _ 1 

20 a2 N 

a 

Les r6sultats pr6c6dents montrent que la variance de la constante est 4 fois plus importante dans I'ap- 
proximation Iin6aire par une constante et 9 fois plus 6lev6e pour I'approximation quadratique et qu*il faut tra- 
vailler sur davantage d'^chantillons d mesure que le nombre de coefficients de I'approxinnatton augmente pour 
conserver les nnemes performances. 

Si les phases corrig^es sont de la forme 

On = A3 n2 + 83 n + C3 

Le choix de la longueur de la fenetre L2 est determine en fonction de I'erreur maximale ef|„ obtenue par 
approximation lin^aire, suivant la relation 

etm = A3{L2 - 1)(L2 - 2)/ 6 
^ Avec une rampe Doppler de dAf Hz par seconde et une dur6e de symbole T, A3 est 6gale d 

A3 = T2 dAf/2 

Dans ces conditions pour obtenir une erreur effn infdrieure d X fois I'^cart type des phases du au bruit a^, il 

suff it pour L2 assez grand que L2 v6rif le la relation 

L ^1 V 12 ^ X" 

^ ^ T dAf V2S/B 

Avec S/B = 3 dB, un rythme symbole de 1000 Hz et une rampe Doppler de 100 Hz par seconde, la relation 
pr§c6dente montre que L2 doit §tre inf^rieur ^ 75 pour que Terreur de iin^arit^ ne d^passe pas le 1/1 Oe de 
I'erreur due au bruit 

Bien que tr6s performant, le precede d^crit pr6cedemment peut encore etre am6llor6. 
^ Parexemple, si d6s la phase de recherche de la position de depart optimale. Terreur quadratique minimale, 

est trap 6lev6e, un autre essai peut encore etre tent6 avec une fenetre de largeur L reduite en esp^rant trouver 
cette fois-ci un r^sultat convenable. 

Comme il subsiste toujours une ambiguity de phase de n dans la determination des phases de reference, 
il n'est pas possible de savoir a priori si la phase du premier echantillon de signal est 0 ou tc puisque le bit 
^ modulant peut avoir Tune ou I'autre valeur 0 ou 1. Cette ambiguity peut dtre levde par tout precede connu, 
comme par exempte celui consistant d insurer une sequence de K bits connue. ou un codage diff6rentiel des 
donnSes avant Emission. 

Bien entendu, le precede et le dispositif qui viennent d'etre d^crits dans te cadre d'une modulation biphase 
s'appliquant aussi bien pour d*autres types de modulation de phase. 
^ Par exemple. pour une modulation h 4 ^tats de phase. appel6e QPSK dans la littSrature, les dispositifs 

1 , 3 et 4 peuvent etre adapt6s en corrigeant les phases et differences de phase dans le dispositif 3 pour les 
ramener entre -«/4 et +fc/4; et pour corriger 6ventuellement Op dans les dispositifs 3 et 4 par pas de 7i/2 et non 
pas n. 

Pour une modulation MSK ou OQPSK. la seule modification d apporter consiste alors d multiplier par j 
^ tous les 6chantillons d'ordre impair par exemple. 

Enf in, le proc^d^ peut aussi §tre adapts pour aff iner la synchronisation du r^cepteur. 

En disposant de K ^chantillons par symbole, au lieu de 1 comme il a et6 suppose jusqu'ici, et done de N.K 
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6chantillons au total, il suffit alors d'effectuer une demodulation complete d partir des 6chantlllons 1 , K+1, ...» 
etc., fournlssant une valeur Emin(1); une demodulation complete ^ partir des 6chantillons 2, K+2, etc.. four- 

nissant une valeur Em,n(2); etc. et une demodulation complete e partir des echantillons K. 2 K etc., four- 

nissant una valeur Emin(K). 

Dans ce cas seul le resultat de la demodulation de rang k qui fournit la plus faible valeur E„in(k) est e retenir 
et le rang k designe alors la position qui est la plus proche de la position de synchronisation ideaie. II est aussi 
possible eventuellement d'effectuer une interpolation de la valeur E^m au voisinage de la k^"^ position pour 
aff iner ce resultat. 



Revendications 



1. Procede de demodulation coherente par bloc de suites de signaux complexes Sn modules en phase sui- 
vant des etats de phase determines caracterise en ce qu'il consiste e extraire (5) la phase de chacun des 

IS signaux complexes en supprimant les sauts de phase inconnus, e rechercher (6) dans la suite des phases 

obtenues une premiere plage de phase de longueur detenminee L pour laquelle la variation de phase est 
la plus lineaire possible, d determiner (7) une deuxieme plage de phases de longueur L2 superieure d L 
correspondant e une approximation lineaire plus precise, e effectuer (8) iterativement des regressions 
lineaires d Taide d'une fenetre gllssante de longueur egale ^ la deuxieme plage de phase L2 pour suppri- 

20 mer les sauts de phase residuels pour effectuer (9) une regression quadratique sur les phases restantes 

et pour calculer une estimation globale des phases de reference des signaux de cheque bloc, et e effec- 
tuer (10) une comparaison de Terreur quadratique moyenne entre les phases de reference estimees et 
les phases brutes obtenues par regression lineaire. 

25 2. Procede selon la revendication 1 caracterise en ce que la suppression (5) des phases inconnues consiste 
pour une modutatk>n e deux etats de phase (O-tt) k corriger les phases des echantillons Sp en les ramenant 
dans rintervalle de phase 



30 
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e calculer les differences de phase successives en les ramenant elles aussi dans rintervalle de phase 



et e reconstituer une suite On de valeurs de phases representatives de la phase de reference du signal. 

3. Procede selon Tune quelconque des revendications 1 et 2 caracterise en ce qu'il consiste pour rechercher 
la premiere plage de phase, e calculer par regression lineaire sur la longueur L I'erreur quadratique moyen- 
ne de phase en prenant successivement chacun des signaux Sn comma origine et e choisir pour origine 
celui qui donne Tenreur quadratique moyenne minimum. 

45 

4. Dispositif pour la mise en oeuvre du precede selon Tune quelconque des revendications 1 e 3 caracterise 
en ce qu'il comporte un dispositif de memorisation (2) des signaux complexes couple e un dispositif (3) 
d'estimation de la phase de reference du signal et un dispositif (4) couple au dispositif (3) d'estimation 
de phase pour reconstituer les echantillons complexes et decoder I'information transmlse par les echan- 

^ lillons. 
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